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Grand- Canonical simulation of 4D simplicial quantum gravity 
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A thorough numerical examination for the field theory of 4D quantum gravity (QG) with a special emphasis on 
the conformal mode dependence has been studied. More clearly than before, we obtain the string susceptibility 
exponent of the partition function by using the Grand-Canonical Monte-Carlo method. Taking thorough care of 
the update method, the simulation is made for 4D Euclidean simplicial manifold coupled to Nx scalar fields and 
Na U(l) gauge fields. The numerical results suggest that 4D simplicial quantum gravity (SQG) can be reached 
to the continuum theory of 4D QG. We discuss the significant property of 4D SQG. 



1. Introduction 

Until now 4D SQG has been investigated from 
several points of view, such as the phase struc- 
ture, nature of the transition, effects of additional 
matter fields, and the modified measure action 
(l). Unfortunately the numerical results could 
not resolve the problem for the existence of the 
continuum theory of 4D QG clearly. Gradually 
the original enthusiasm to explore 4D QG has 
been toned down. However, we remind that many 
problems about 4D SQG are yet to be studied. 

As one of those problems, we pay attention to 
additional matter effects in 4D Euclidean SQG 
(Eucl.SQG). According to the recent numerical 
simulations, the string susceptibility exponent 
( 7 ( 4 )) takes a negative value with adding a few 
U(l) gauge fields. It makes drastic change of the 
phase diagram, and continuous phase transition 
with one gauge matter has been found ||. More- 
over, the numerically obtained 7^ 4 -* as a function 
of the number of matter fields turned out to be 
approximately equal to that of the prediction of 
4D conformal field theory ||. 

In order to see more definitely how the matter 
fields act on 4D Eucl.SQG, we repeat the study 
with employing the Grand-Canonical method 



used for 2D SQG[|, which has given the satis- 
factory agreement of two results. For discussing 
the continuum theory of 4D QG based on the nu- 
merical analysis, we pay attention to the func- 
tional form of the partition function, which we 
expect enough to gain information for selecting 
the correct theory of gravity, if it exists at all. 
For obtaining the functional form of the partition 
function, we need to perform a simulation with 
varying the space volume. In practice, we calcu- 
late the probability (P(A^4)) as a function of the 
volume (N4) by the Grand-Canonical method of 
the dynamical triangulation (DT) and estimate 
7" from the scaling behavior of P(A r 4). 

2. Model and Grand-Canonical method 

Let us explain our model and the numerical 
method of the Grand-Canonical simulation. We 
consider a system of the simplicial gravity coupled 
to U(l) gauge fields and massless scalar fields. 
The action (5) is given on 4D simplicial manifold 
as 



S = Sq + Sa + Sx , 



(1) 



* presented by S.Horata 



where Sq, Sa and Sx denote the action for the 
gravity U(l) gauge fields (A) and scalar fields 
(X), respectively. 



For the gravity part, we use the discretized 
Einstein-Hilbert action in 4D, 



S g [k 2 , k 4 ] = k 4 N a - n 2 N 2 , 



(2) 



where Ni denotes the number of i-simplex. The 
two parameters, k 2 and K4, correspond to the in- 
verse of the gravitational constant and the cosmo- 
logical constant, respectively. The action for the 
Nx scalar fields X", a — l...Nx, on the vertex i 
is given by 



(3) 



o=l ij 



where o{Uj ) is the number of four-simplices shar- 
ing a link lij. And the action for the Na vector 
field Afj( a — 1...Na) on link 1^ reads 



a=l ijk 

where o(tyfe) is the number of four-simplices shar- 
ing a triangle fcyfc . These actions give the correct 
thermo-dynamical limit for large N4 because they 
are proportional to N4. 

In order to simplify the calculation, the Monte- 
Carlo simulations have been performed with the 
Canonical method previously. This method re- 
quires that N4 fluctuates around the given target 
volume (V4). Thus, the volume fixing term (SS) 
is added to the total action Eq. ([!]), 



SSiNi, V4) = <5k 4 - V 4 y , 



(5) 



where Sk 4 is the free parameter. However, the 
Canonical method has limitations about the finite 
size effect and the ergodic property. Furthermore, 
the functional form of the size dependence of the 
partition function cannot be given clearly. For the 
improvement of numerical accuracy we perform 
the simulations varying the space volume by the 
Grand-Canonical method. The Grand-Canonical 
partition function with the matter fields is given 
as 



Z(k 2 ,k 4 ,N a ,N x ) = ^2 1 



n /n^ 

N A ,N X V l£T 
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-Sg(«2i«4) 



(6) 



where all geometries T and matter fields A, X are 
updated using the Metropolis method keeping the 
topology of discretized surface to be S 4 . 

From the analogy of the 2D case , we assume 
that the partition function obeys the KPZ like 
formula for the large- N4 limit as 



7 W («2)-3 e ( K4 - re J(K 2 ))JV 4> 



(7) 



N.i 



Then the probability P(A^) for the large-7V4 be- 
haves approximately as 



P(Ni) oc iY2 <4) -V (K4 -^ )Ar4 



(8) 



By tuning K4 close to k 4 , the exponential growth 
is suppressed and F(A r 4 ) is given as the power 
function of N±. The parameter K4 is determined 
in the following manner: If K4 is chosen to be 
less (or greater) than the system will expand 
(or shrink) exponentially in the Grand-Canonical 
simulation. By requiring an exponential increase 
(or decrease) to disappear in the simulation, we 
can find K4 precisely with the following form, 

After tuning K4, we generate the configurations 
constraining N4 within an upper- and a lower- 
bound (7V 4 min < iV 4 < N r 4 nax ). 

3. Numerical result 

In this section, we report numerical results with 
the Grand-Canonical method. In order to ob- 
tain P(Ni), we allow N4 to fluctuate between 

N min = 1Q00 and N ™ax = 10000 j n Fig g we 

show P(N 4 ) versus N4 at K4 close to k\. The 
numerical result shows that P{N4) is given as a 
power function of N4. It is suggested that the 
partition function in 4D obeys the KPZ type par- 
tition function as the 2D case. 

We can also compute 7 < - 4 - ) from the slope of 
P(N4). Then, we show the parameter (6), which 
is calculated from j 4 , 




(10) 
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Figure 1. We plot P(N^) versus N 4 with log-log 
scale. 

This relation is given from the analogy of 2D 
QG and the parameter b is similar to the cen- 
tral charge in 2D QG. In Fig]|, we plot b ver- 
sus the number of matter fields, Nx + Q2Na, 
on the critical point (n^) between the crumpled 
phase and the smooth phase. We compare b ob- 
tained by two different numerical methods, i.e. 
the Canonical method with the MINBU algo- 
rithm and the Grand-Canonical method. Error 
bars include both the errors in determining n\ 
and the statistical errors, and we found that the 
former method has much larger errors than the 
latter method. We made the parameter fit of the 
data of the Grand-Canonical method by a linear 
function, b = 0.0030(3) {N x + 62N A ) + 3.98(3). 
The slope b is approximately equal to the analyt- 
ical result! 



4. Summary and Discussion 

Let us summarize and discuss the main points 
made in the previous section. From the numerical 
results of the Grand-Canonical method, we find 
the KPZ type partition function also appears in 
4D QG. The string susceptibility exponent 7^ 
obtained from the size dependence of the par- 
tition function can be obtained by the Grand- 
Canonical method rather accurately. The matter 
contribution to 6, which is computed from 7" 
with the analogy of 2D QG, turns out to be ap- 
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Figure 2. We plot b versus the number of the 
matter fields, N x + 62N A - 



proximately equal to the conformal part [|5j|6j] . Our 
numerical results suggest that it may be possible 
to understand the continuum theory of 4D QG 
from the analogy of 2D QG. In order to look more 
details on the relation to the continuum theory of 
4D QG, these results encourage us to continue in- 
vestigation by the Grand-Canonical method and 
push more efforts on the analytical calculation of 
the conformal gravity at the same time. 
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